Draft version February 2, 2008 

Preprint typeset using L^T^iX style emulateapj v. 1 1/12/01 



CONSTRAINTS ON EXTENDED QUINTESSENCE FROM HIGH-REDSHIFT SUPERNOVAE 
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ABSTRACT 

We obtain constraints on quintessence models from magnitude-redshift measurements of 176 type la Super- 
novae. The considered quintessence models are ordinary quintessence, with Ratra-Peebles and SUGRA poten- 
tials, and extended quintessence with a Ratra-Peebles potential. We compute confidence regions in the il m Q-a 
plane and find that for SUGRA potentials it is not possible to obtain useful constraints on these parameters; for 
the Ratra-Peebles case, both for the extended and ordinary quintessence we find a <L0.8, at the la level. We also 
consider simulated dataset for the SNAP satellite for the same models: again, for a SUGRA potential it will not 
be possible to obtain constraints on a, while with a Ratra-Peebles potential its value will be determined with an 
error ^,0.6. We evaluate the inaccuracy made by approximating the time evolution of the equation of state with 
a linear or constant w(z), instead of using its exact redshift evolution. Finally we discuss the effects of different 
systematic errors in the determination of quintessence parameters. 

Subject headings: Cosmology: observations, theory, cosmological parameters - Supernovae: general 



1. INTRODUCTION 

Cosmological tests such as the baryon fraction in galaxy clus- 
ters (Hradecky et al. 2000; Allen et al. 2002), the abundance of 
massive galaxy clusters (see e.g. Bahcall 2000; Borgani et al. 
2001), the magnitude-redshift relation for type la Supernovae 
(Riess et al. 1998; Perlmutter et al. 1999; Gott et al. 2001; 
Tonry et al. 2003), the statistical analysis of the galaxy dis- 
tribution in large redshift catalogs (e.g. Percival et al. 2001; 
Verde et al. 2002) and Cosmic Microwave Background (CMB) 
anisotropies (Bennett et al. 2003; Spergel et al. 2003) indicate 
a low value for the matter density parameter today, fi m o, proba- 
bly lying in the interval 0.15 <^£l„,o ^0.4. Recent results from 
the study of CMB anisotropies obtained with the WMAP satel- 
lite (Bennett et al. 2003) also provide strong support for a flat 
(or very nearly flat) Universe. Combining these two different 
indications leads to the hypothesis that a new form of energy, 
named dark energy (DE), fills the gap between f2 m o and unity: 

fl m + floEO = i ■ 

One of the main goals of modern cosmology is to explain the 
nature of dark energy. The simplest solution to this problem is 
to introduce a cosmological constant A in our Universe model. 
This scenario can have a simple theoretical interpretation, as A 
can be related to the energy density of the vacuum state in quan- 
tum field theory (Zel'dovich 1968). Unfortunately, this simple 
explanation results in a very large (formally infinite) value for 
the vacuum energy density, which is larger by tens of orders of 
magnitude than the observed one (of the order of 10 -47 GeV 4 ). 
What emerges is a fine-tuning problem, namely the "cosmolog- 
ical constant problem". Another apparently unnatural feature 
of the cosmological constant model is that it starts to dominate 
the Universe evolution only in the very near past. This issue 
is usually referred to as the "cosmic coincidence problem" and 
reduces to a fine-tuning problem in the choice of the initial con- 
ditions, in particular of the value of A. 

A possible way of alleviating these problems is to allow for 
a time variation of the dark energy density, which is constant 
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if due to A. A very interesting class of models with this prop- 
erty goes under the name of quintessence (Coble et al. 1997; 
Caldwell et al. 1998; Ferreira & Joyce 1998; Viana & Liddle 
1998). 

One of the most promising cosmological tests on the prop- 
erties of the dark energy component is based on the already 
mentioned magnitude-redshift relation for type la Supernovae 
(see e.g. Brax & Martin 1999; Podariu & Ratra 2000; Podariu 
et al. 2001; Goliath et al. 2001; Weller & Albrecht 2001, 2002; 
Eriksson & Amanullah 2002; Gerke & Efstathiou 2002; Pad- 
manabhan & Choudhury 2003; Di Pietro & Claeskens 2003; 
Knop et al. 2003). In fact, these objects can be considered as 
good standard candles, which makes it possible to determine 
their luminosity distance, whose dependence on redshift is spe- 
cific of each particular cosmological framework. 

In this paper we will focus on three different kinds of quin- 
tessence models whose features are briefly described in Sec- 
tion 2. In Section 3 we present the constraints obtained on 
these models using the magnitude-redshift relation for existing 
type la Supernovae (SNIa) data. In Section 4 we illustrate how 
the constraints will improve with the SuperNovae Acceleration 
Probe (SNAP) satellite which is currently being projected and 
will be devoted to the discovery and study of SNIa (Aldering et 
al. 2002; see also: http://snap.lbl.gov); we also check the va- 
lidity of approximating the exact redshift evolution of the equa- 
tion of state with a linear or constant behavior and discuss the 
possible effects of different systematic errors on the parameter 
determination. Finally, in Section 5 we present our conclusions. 

2. THEORETICAL FRAMEWORK 

In this paper we focus on the Extended Quintessence (EQ) 
models, introduced in (Perrotta et al. 1999; Baccigalupi et 
al. 2000). Extended quintessence and related models have 
been also considered in (Chiba 1999, 2001; Uzan 1999; Bar- 
tolo & Pietroni 2000; Boisseau et al. 2000; de Ritis et al. 
2000; Faraoni 2000; Fujii 2000; Chen et al. 2001; Bean 2001; 
Gasperini 2001 ; Perrotta & Baccigalupi 2002; Riazuelo & Uzan 
2002; Torres 2002; Kneller & Steigman 2003). 

For these models the evolution of the scale factor a and the 
scalar field (j> responsible for the quintessence component can 
be obtained by solving the set of equations: 
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In the previous equations the dots denote derivatives with re- 
spect to the conformal time and the subscript denotes differ- 
entiation with respect to the scalar field; R is the Ricci scalar; 
P fluid represents the energy density associated with all the con- 
stituents of the Universe except for the quintessence scalar field; 
V ((/)) is the quintessence potential and finally F (</>) is a func- 
tion specifying the form of the coupling between (f> and gravity. 
Hereafter we will always refer to the non-minimally coupled 
(NMC) scalar field models (Perrotta et al. 1999), for which the 
function F is defined as F{4>) = 1 /8wG+ F {4>) — F ((/>o), where 
F((j>) = £,(j) 2 . This kind of models has two free parameters: 
the dimensionless constant £, parametrizing the amount of cou- 
pling, and the present value of the scalar field <po- 

The coupling between the scalar field and gravity generates 
a time-varying gravitational constant (see e.g. the review by 
Uzan 2003). Upper bounds on this variation come from lo- 
cal laboratory and solar system experiments (Gillies 1997) and 
from the effects induced on photon trajectories (Reasenberg et 
al. 1979; Will 1984; Damour 1998). As pointed out by Perrotta 
et al. (2000), these bounds become constraints on the parame- 
ters of the models: 
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We will use the previous inequality in the next sections in 
order to improve our determination of the cosmological param- 
eters. 

For completeness and in order to allow a comparison with 
similar analyses, we will also consider the case of ordinary 
quintessence (OQ), i.e. models for which there is no direct cou- 
pling between the scalar field and gravity (it is often referred to 
as the minimal coupling case). OQ can be easily obtained from 
EQ in the limit of £ -> 0. 

If we want to completely specify a quintessence model, we 
have to choose the analytical form for the potential V ((/)). One 
of the main advantages of a time-varying dark energy density 
is that it is possible to alleviate the cosmic coincidence prob- 
lem. This is achieved by assuming particular classes of poten- 
tials, the so-called "tracker potentials" (Steinhardt et al. 1999), 
for which one obtains at low redshifts the same time evolution 
for a and <j>, even starting from initial conditions which differ 
by orders of magnitude: this leads to a dark energy dominated 
era close to the present time, without fine-tuning on the initial 
conditions. For the following analysis, we will consider two 
different classes of tracker potentials: the inverse power-law 
Ratra-Peebles potential (hereafter RP; Ratra & Peebles 1988; 
see also Peebles & Ratra 2003 and references therein) 

vw = -jsr, (4) 

and the SUGRA potential (Binetruy 1999; Brax & Martin 1999): 

M A+a 

y = — exp(47rG0 2 ). (5) 

In particular, we will use the potential (4) in the context of both 
EQ and OQ models, while the potential (5) will be considered 
in the minimal coupling (£ = 0) case, only. 

We solved numerically Eqs.(l) and (2) in the tracking regime. 
The behaviors we found for a(t) and <p(t) (not shown in figure) 
are in excellent agreement with those obtained by Baccigalupi 
et al. (2000), whose analysis also includes the Ratra-Peebles 
case in the framework of EQ. In the next section we will use 
these results to obtain our theoretical estimates of the luminos- 
ity distance di- 



3. CONSTRAINTS FROM PRESENT HIGH-REDSHIFT 
SUPERNOVAE DATA 

The purpose of this section is to test the possibility of con- 
straining the cosmological parameters describing quintessence 
models by using the best SNIa dataset presently available. To 
this aim we build a sample combining data coming from the 
literature. As a starting point, we consider the data reported 
in Table 15 of Tonry et al. (2003). In particular, we use for 
our analysis the subset which is presented by the authors as the 
most reliable one for cosmological studies. This data compi- 
lation, comprising 172 SNIa, is obtained by eliminating from 
the original whole sample of 230 SNIa the objects at low red- 
shift (z < 0.01) and those with high reddening (Ay > 0.5 mag). 
Then, we consider the data from Table 3 of Knop et al. (2003), 
but including in the sample only SNIa belonging to their "low- 
extinction primary subset" (7 objects). In their Table 4, Knop 
et al. (2003) present the new fits to the Perlmutter et al. (1999) 
data, which are already included in the Tonry et al. (2003) 
sample. In this case we decided to use the magnitudes from 
Knop et al. (2003), because they are obtained using new fitting 
lightcurves. The two catalogues have also SNIa in common in 
the low-redshift sample, taken from Hamuy et al. (1996) and 
Riess et al. (1999). For these we prefer to use the data re- 
ported by Tonry et al. (2003) because they are obtained as a 
median of different fitting methods. For coherence with our 
previous choice, we decided to exclude 5 objects, present in the 
low-redshift and Perlmutter et al. (1999) samples, but excluded 
by the "low-extinction primary subset" of Knop et al. (2003). 
Finally, we add two new SNIa, 2002dc and 2002dd, recently 
studied by Blakeslee et al. (2003). Therefore the sample of 
SNIa here considered comprises 176 objects. 

To constrain the cosmological parameters, we compare through 
a x 2 analysis the redshift dependence of the observational es- 
timates of log (di) to their theoretical values, which for a flat, 
matter or dark energy dominated universe can be obtained as 

-,-1/2 

ft m0 (i+O +n 4>0 f( Z , y 
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where 



f(z) = exp ( 3 / — dz 



dz', (6) 
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Here Ho is the present value of the Hubble constant, and O m o 
and O^o represent the contributions to the present density pa- 
rameter due to matter and scalar field, respectively. The quan- 
tity w = p^jp$ sets the equation of state relating the scalar field 
pressure 

6 2 3HF 

r* = h +v W— <*> 

and its energy density 

12 p ffp 

(9) 
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The quantity w(z) can be computed by using the numerical so- 
lutions of the previous equations, once a cosmological model is 
assumed. 

We can now find, for the different quintessence models de- 
scribed in the previous section, the parameters (and their con- 
fidence regions) which best fit the SNIa data by using the stan- 
dard x 2 method. For this analysis we have to consider the er- 
rors on the distance moduli. For the objects coming from Tonry 
et al. (2003), they are obtained directly from their Table 15, 



3 




0.2 0.4 0.2 0.4 0.6 



^m0 ^mO 

FIG. 1. — Confidence regions (A\ 2 = 2.30,5.99, 11.8, corresponding to 1, 
2 and 3cr for a Gaussian distribution with two free parameters, respectively) 
for the parameters f!,„o and a, as obtained using the sample of 176 SNIa 
data presently available. Different panels refer to different quintessence mod- 
els: ordinary quintessence with a Ratra-Peebles potential (top left); ordinary 
quintessence with a SUGRA potential (top right); extended quintessence with 
a Ratra-Peebles potential (bottom left); extended quintessence with a Ratra- 
Peebles potential when the upper limits on the time variation of the gravita- 
tional constant are satisfied (bottom right). 



once a value for Hq is assumed. On the other hand, Knop et al. 
(2003) report errors for the apparent magnitudes of their SNIa. 
We obtain the uncertainty on the distance modulus by adding 
in quadrature an error of 0.05 in the estimate of the absolute 
magnitude, as suggested by the data of Hamuy et al. (1996), 
Riess et al. (1999) and Knop et al. (2003), Table 4. In addi- 
tion, we include an extra contribution to the error coming from 
the possible uncertainty in the peculiar velocities. In particular, 
following Tonry et al. (2003) we add 500 km s" 1 divided by the 
redshift in quadrature to the distance error. 

We compute the values for x 2 on a regular grid of the con- 
sidered parameters (Q m o, a and £ for the EQ models; fi mu and 
a for the OQ models), looking for its minimum value, xLin- ^ n 
particular, we allow fl m o to vary between and 1 with spacing 
of 0.001; a between and 12 with spacing of 0.1 and £ between 
0.001 to 0.100 with spacing of 0.001. The results in the plane 
ttmQ-ct are shown in Figure 1, where we display the contours of 
constant A\ 2 = X 2 ~X 2 B m = 2.30,5.99 and 11.80, correspond- 
ing to lcr, 2er and 3er for a Gaussian distribution with two free 
parameters, respectively. 

The top-left panel refers to the ordinary quintessence mo- 
del with a RP potential. In this case we find Xmin = 206 for 
^mO = 0.30 and a = 0.0. However, from the plot it is evident 



that there is a strong degeneracy between the two free parame- 
ters a and fi mu and it is not possible to obtain strong constraints 
on them at the same time. Nevertheless, we can extract some 
information by considering the x 2 distribution when only a sin- 
gle parameter is allowed to vary. In this case, we obtain the 
errorbars associated to the best fit value by assuming Ax 2 = 1, 
which corresponds to Id for a Gaussian distribution with one 
single parameter. The SNIa dataset we used does not allow to 
obtain tight constraints on the parameter a: a < 0.83 at the lcr 
confidence level, with a best fit value of a = 0. In addition, we 
obtain il„,o = 0.30^ jq, for the matter density parameter. If we 
impose the Gaussian prior fi m o = 0.27 ± 0.04, as suggested by 
the combined analysis of recent CMB observations and large- 
scale structure data (Spergel et al. 2003), we obtain a < 0.47 at 
the lcr confidence level (best fit value: a = 0. 10). We notice that 
our results are in agreement with those obtained from a similar 
analysis carried out by Podariu & Ratra (2000), even though, 
thanks to the improvement in the SNIa dataset, the confidence 
contours start to close off, at least at the lcr confidence level. 

In the top-right panel we report the results still for the OQ 
model, but with the SUGRA potential. The values for Xmm an d 
the corresponding parameters are the same obtained for the RP 
potential. In this case the confidence regions are almost vertical, 
showing a very small dependence on a: this does not allow us 
to extract constraints on fi m o an d a - Only if we impose the 
Gaussian prior £l„,o = 0.27 ± 0.04, we are able to obtain a < 
2.78 at the la level (best fit value: a = 0.24). 

In the two bottom panels of Figure 1 we present the results 
for the model with EQ and a RP potential. In this case we 
have three free parameters: in addition to a and il,„o, there is 
£, which parametrizes the strength of the coupling between the 
scalar field and gravity. In the bottom-left panel we show the 
two-dimensional confidence regions in the a - £!„,o plane as ob- 
tained by minimizing \ 2 with respect to £: they appear very 
similar to the previous case of ordinary quintessence with RP 
potential, only slightly larger. Again it is convenient to discuss 
the results when a single free parameter is considered. Un- 
fortunately, there are no possibilities to obtain constraints on 
£ (best fit for £ = 0.001), even at the lcr confidence level. For 
the other two parameters our results are again very similar to 
that of ordinary quintessence with RP potential: a < 0.82 (best 
fit a = 0.02) and fi m o = 0.28+oQg. The bottom-right panel shows 
how the confidence regions change when the upper limits on 
the time variation of the gravitational constant are taken into 
account. The constraints on a and f2 mu do not change signifi- 
cantly: a < 0.82 and f2 mu = 0.28^Q gg, respectively. Considering 
£, the inclusion of the upper bounds does not improve the situ- 
ation, and £ remains completely undetermined. 

Finally, instead of imposing the constraints just described, 
we can consider the prior £l,„o = 0.27 ± 0.04. Again, we cannot 
obtain useful bounds on £, while we obtain a < 0.45 at the lcr 
level (best fit: a = 0.1 1). The previous results clearly show the 
necessity for new samples of SNIa. The SNAP satellite will 
be dedicated to this purpose: in the next section we will check 
the improvement that it will allow on the determination of the 
cosmological parameters. 

The results presented in this section assume that all the data 
points are purely governed by statistical errors. Even if we 
carefully selected our SNIa sample in order to reduce possible 
systematic errors, their importance in the present data cannot 
be completely excluded. As largely discussed in Knop et al. 
(2003, see also Tonry et al. 2003), systematics can have dif- 
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ferent origins, In particular Knop et al. (2003) studied the ef- 
fects on the determination of the cosmological parameters (they 
considered fi m rj, Q^o, w) coming from different SNIa selection, 
lightcurve fitting methods, contamination from non-type la su- 
pernovae, Malmquist bias, gravitational lensing, dust proper- 
ties, etc. Their conclusion is that the total systematic error is of 
the same order of magnitude as the statistical uncertainty. Be- 
cause of the weakness of the constraints we obtained on a and £ 
using present SNIa data, we prefer to leave an extended discus- 
sion of systematics to the next section, where we will present 
the expected results from the SNAP satellite. However, here 
we report the results of some tests which have been done. In 
order to verify if our sample suffers from selection effects, we 
repeated our analysis by considering only 54 SNIa belonging 
to the low-extinction primary subset of Knop et al. (2003) and 
coming only from the Supernova Cosmology Project. Because 
of the smallness of this sample, the resulting constraints cannot 
be directly compared to those obtained from the whole sample 
of 176 SNIa, being weaker, but we notice that the best-fit val- 
ues for O m o and a are only slightly different. Then, we checked 
the systematic effect due to the possible type-contamination, 
by excluding from our sample those objects whose confirma- 
tions as type la supernovae are questionable. This new sub- 
sample has a smaller number of very-high redshift SNIa, and 
consequently the constraint on a becomes much weaker when 
compared with the whole sample: a < 1.62 for OQ with RP 
potential, and a < 1 .52 for EQ, always with RP potential. Fi- 
nally, we checked that the exclusion of very-low redshift objects 
(z<0.03), whose measurements of distance moduli are possibly 
affected by peculiar velocities, does not change the resulting 
confidence levels, confirming that our treatment of this kind of 
error is reliable and that the constraining power of the analysis 
comes from high-redshift SNIa. 

4. CONSTRAINTS FROM FUTURE HIGH-REDSHIFT 
SUPERNOVAE DATA 

In order to improve the sample of studied SNIa, a new satel- 
lite, the SuperNovae Acceleration Probe (SNAP) is currently 
under project. SNAP is expected to perform in two years a 
complete spectroscopic and photometric analysis for approxi- 
mately 2000 high-redshift SNIa reaching a maximum redshift 
of z = 1.7. Due to the large increase of both the number of 
observed objects and the covered redshift interval, SNAP can 
prove extremely useful for cosmological studies. 

In order to check the ability of this satellite to constrain the 
parameters for quintessence models, we generate pairs of dis- 
tance moduli and redshift for 1998 SNIa, adopting the fiducial 
distribution proposed by Kim et al. (2004) and shown in Figure 
2. Always as in Kim et al. (2004), we include in addition a 
very low-redshift sample of 300 supernovae with redshift uni- 
formly distributed between z = 0.03 and z = 0.08. Thus, the total 
number of objects considered in the following analysis is 2298. 

The simulated distance moduli fi are assumed to have a Gaus- 
sian distribution around the true value, with a dispersion ct m = 
0.15. The true distance modulus for a given redshift z is known 
once we specify the background cosmology. For our analysis, 
we decided to simulate SNIa datasets with three different back- 
ground cosmologies, whose properties and parameter values are 
reported in Table 1 . For all datasets, we use for the Hubble con- 
stant Ho = 71 km s" 1 Mpc" 1 , which is the value suggested by 
the combined analysis of WMAP and large-scale structure data 
(Spergel et al. 2003). 
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FIG. 2. — The SNIa redshift distribution used in the simulated SNAP samples 
(see also Kim et al. 2004). 



Table 1 

Properties and values of the parameters for the 
different cosmological models assumed in the 
generation of simulated snap samples. 



Name 


Model 


Potential 


a 


^m0 




dataset 1 


OQ 


RP 


1.0 


0.3 


0.000 


dataset 2 


OQ 


SUGRA 


1.0 


0.3 


0.000 


dataset 3 


EQ 


RP 


1.0 


0.3 


0.015 



Adopting the same x 2 method described in the previous sec- 
tion, we then analyze these simulated datasets, fitting each quin- 
tessence model only with the data obtained from the same back- 
ground cosmology. In computing the value for x 2 we consider 
an error of 0.15 mag on the distance moduli and we neglect the 
errors on the redshifts z- 

It is also quite important to verify how good is the approxi- 
mation of using a linear (i.e. w(z) = wo + wiz, with wq and w\ 
suitable constants) or a constant (w = Wo = const.) equation of 
state in the fitting procedure (as often done in the literature), 
rather than its exact redshift evolution (as we did in the pre- 
vious section). To this purpose we apply the following proce- 
dure. Let us refer, for simplicity, to one of the OQ models. 
For each pair of values (O m o, a), we numerically obtain the red- 
shift evolution of the equation of state computing w(z) in 1700 
equally spaced values of z in the interval [0, 1.7]. Using the 
Least squares method, we determine the straight line which best 
fits this ensemble of points. Finally, we use the theoretical dis- 
tance modulus obtained from this linear approximation for w 
instead of its exact evolution to calculate the value of x 2 with 
the SNIa simulated data. Similarly, in order to check the va- 
lidity of the approximation of a constant equation of state, we 
follow the same procedure but substituting the straight line with 
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FIG. 3. — Confidence regions (Ax 2 = 2.30,5.99, 11.8, corresponding to 1, 2 
and 3<t for a Gaussian distribution with two free parameters, respectively) for 
the parameters Q,„q and a for the ordinary quintessence model with a Ratra- 
Peebles potential. The results are obtained with the simulated SNAP sample 
(dataset 1). Panel (a) refers to the constraints obtained by assuming the exact 
equation of state »'(z); panel (b) to the linear w(z) approximation; panel (c) to 
the constant w approximation. The final panel shows the superposition of the 
three cases. 



the mean value of w(z) in the considered range. In the follow- 
ing figures we will illustrate the differences in the confidence 
regions between these two approximations and the results ob- 
tained assuming the exact equation of state. 

In Figure 3 we show the confidence regions in the a - ft m o 
plane obtained by fitting the OQ model with a RP potential to 
the dataset 1 . As done also for the other fits described in this 
section, we fix the value of the Hubble constant to the 'true' 
value, Ho = 71 km s" 1 Mpc" 1 . In the panel (a) we show the re- 
sults obtained with the exact equation of state. The best-fitting 
parameters are: a = 0.91^33 and n ml) = 0.30^;^ at the la level. 
They are consistent with the values assumed in the background 
cosmology, within the error limits. By comparing these results 
to those presented in the top-left panel of Figure 1 , it is also ev- 
ident the improvement with respect to the constraints obtained 
with the present SNIa data: the confidence regions are much 
narrower. Nonetheless, there is still a little degeneracy between 
a and £l,„o: to solve it will be useful to combine magnitude- 
redshift measurements to other cosmological tests (see, e.g., 
Balbi et al. 2003; Frieman et al. 2003; Jimenez et al. 2003; 
Caldwell & Doran 2003). The case of OQ model with a RP po- 
tential has already been examined by Podariu et al. (2001) for 
SNAP simulated data. Their confidence regions in the il„,o - a 
plane look quite similar to ours, even if they considered a fidu- 



cial cosmology with a = 4 and fl m o = 0.2. A similar analysis has 
been carried out also by Eriksson & Amanullah (2002), who 
found an error on a which is approximately twice larger than 
ours. This difference is probably due to the different redshift 
distribution adopted. 

In panel (b) of Figure 3 we show the confidence regions ob- 
tained when x 2 is computed using the linear approximation 
for w(z), while in panel (c) we adopt a constant equation of 
state. The differences between the three considered cases (ex- 
act, linear and constant w) are evident also from the last panel, 
showing the superposition of the corresponding confidence re- 
gions. The best-fitting values we obtain are a = 0.8 l^i} an d 
O m0 = 0.3 1^ for the linear w(z), and a = 0.76^ and fi m0 = 
0.31^q{]3 for the constant w, in both cases at la level. The mea- 
sured cosmological parameters are then consistent with the true 
values within the error limits. However, there seems to be a 
systematic tendency to underestimate a, which increases when 
we decrease the accuracy used to describe the redshift evolution 
of w. Nevertheless, this systematic offset will be too small to 
be put in evidence by the SNAP data, even at the la level. On 
the contrary, the approximations on w{z) do not strongly influ- 
ence the determination of fi m o. Finally, we notice that, at the 
lcr confidence level, the size of the errors for both a and f2 m o 
appear to be smaller, i.e. systematically underestimated, when 
the linear and constant approximations are used instead of the 
exact solution. 

We then perform the same kind of analysis by fitting OQ 
models with a SUGRA potential using the dataset 2. The re- 
sults are shown in Figure 4. Looking at panel (a), which refers 
to the results obtained with the exact w(z), it is evident that the 
SUGRA models present a strong degeneracy in a, becoming 
larger with increasing values of a. This was already pointed 
out by Brax & Martin (1999), who studied the dependence of 
wo on a in the same class of models. As a consequence, in this 
case it is impossible to determine the errors on the cosmological 
parameters. The confidence regions under the approximation of 
linear equation of state are shown in panel (b): they appear very 
similar to those displayed in panel (a). Instead, if we assume a 
constant w (see panel c), the lcr contour closes off, and we ob- 
tain: a = 0.34^ 33 and f2„,o = 0.32 ±0.02, which are marginally 
consistent with the true values a = 1 and J7„,o = 0.30. Again, the 
last panel shows the superposition of the three cases. The situ- 
ation slightly improves if we impose a Gaussian prior on f2,„o, 
namely O m0 = 0.30 ± 0.05. In fact we obtain a = 0.95!^;^ for 
the exact w(z), a = 0.49! % for the linear w(z) and a = 0.38!° ;^ 
for the constant w, all at the la confidence level. All these 
values are consistent within the error limits with the original 
assumption of a = 1, but we see that by using the exact w(z) 
the determination is very poor. Moreover, the use of linear or 
constant equation of state approximations could lead to an ar- 
tificially higher precision in the determination of a. For these 
reasons, we can conclude that even with SNAP it will not be 
possible to constrain OQ models with a SUGRA potential. 

Finally, we use the dataset 3 to estimate the best-fitting pa- 
rameters for EQ models with a RP potential. To this aim we 
determine the values of \ 2 on a three-dimensional grid of val- 
ues for the parameters a,fl m o and £. In particular, we consider 
100 values for £, regularly spaced between 0.001 to 0.100, as 
done in the previous section with the presently available SNIa 
data. 

In Figure 5 we show the confidence regions in the J7„,o - a 
plane, obtained by minimizing (for each pair of these parame- 
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FIG. 4. — As Figure 3, but for the ordinary quintessence model with a 
SUGRA potential. The results are obtained with the simulated SNAP sam- 
ple (dataset 2). 
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FIG . 5 . — As Figure 3, but for the extended quintessence model with a Ratra- 
Peebles potential. The results are obtained with the simulated SNAP sample 
(dataset 3). 



ters) x 2 with respect to £. Panel (a) refers to the exact w(z) case 
and can be directly compared to the bottom-left panel of Figure 
1 . Even if the regions become narrower, they are still larger than 
in the OQ models. It is possible to obtain the one-dimensional 
distribution separately for each of the three parameters, by min- 
imizing \ 2 with respect to the others. In this way we find that 
there is very little dependence of the minimized \ 2 on £. Then, 
it is not possible to obtain useful constraints on this parameter 
(best fit value: £ = 0.001). Instead, we have a = 0.95^ and 
n„,o = 0.30±o y3 at the lcr level. This values are consistent with 
the true values (a = 1 and £! m o = 0.3) within the errors. 

The confidence regions obtained using the approximated (lin- 
ear and constant) equations of state (panels (b) and (c), respec- 
tively) are quite different from the previous ones, as it is evident 
from the superposition in the last panel. Again, there is no pos- 
sibility to constrain £ (best fit values: £ = 0.001 for linear w(z), 
and £ = 0.003 for constant w). For the remaining parameters 
we obtain: a = 0.84^;^ and n mQ = 0.31^ for the linear w(z); 
and a = 0.78^ 4g and fl m o = 0.32^qj for constant w, at la con- 
fidence level. These values are always consistent with the true 
ones, within the error limits. 

As a final point, it is interesting to discuss how the confi- 
dence regions change if we impose the constraint coming from 
the time variation of the gravitational constant. The results are 
shown in Figure 6, where we draw the confidence regions in the 
f2 m o - a plane, obtained by minimizing % 2 with respect to £ and 
considering only the combination of the parameters for which 
the bounds on G(f ) are satisfied. This condition excludes com- 



bination of the parameters for which a and £ are high, and this 
is the main motivation for the differences with respect to Figure 
5. Considering the case where the exact equation of state is used 
(panel (a)), we have a = 0.95 ± 0.43 and Q m0 = 0.30 ± 0.03 at 
lcr: there is no significant change in the determination of these 
parameters with the imposition of the constraint. However, this 
time we are able to obtain an upper limit on £: £ < 0.028 at the 
lcr level. This result is consistent with the assumed true value: 
£ = 0.015. Panel (b) refers to the linear w(z) approximation, for 
which we derive a = 0.84^J|, Q m0 = 0.3 1^ and £ < 0.030 at 
the lcr level. Finally, panel (c) considers the constant w approx- 
imation. In this case the resulting constraints are a = 0.78^ 33, 
Q m0 = 0.32!°°' and £ < 0.030 at the lcr level. Again, all this 
values are consistent with the true ones. 

Up to now, we discussed the possibility of determining the 
quintessence parameters in the future data coming from SNAP 
by assuming intrinsic statistical errors only. Here we consider 
the impact of systematic errors on the previous results. Kim 
et al. (2004) define two general forms for them: uncorrelated 
systematic uncertainties and magnitude offsets. The first case 
represents a random dispersion which cannot be reduced below 
a given magnitude error over a finite redshift bin (here assumed 
to be Az = 0.1): possible examples are calibration errors and 
imperfect galaxy subtraction. The second case is a coherent 
shift acting as a bias on all SNIa magnitudes: selection effects 
as Malmquist bias or detector problems can produce this kind 
of effect. 

In order to simulate irreducible systematics, we strictly fol- 
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FIG. 6. — As Figure 3, but for the extended quintessence model with a 
Ratra-Peebles potential. The results are obtained by using the simulated SNAP 
sample (dataset 3) and imposing the upper bound on the time variation of the 
gravitational constant. 



low Kim et al. (2004). In particular, for each redshift bin, 
containing Nsn objects (see Figure 2), we add in quadrature 
to the intrinsic magnitude dispersion per SNIa (0.15 mag) an 
irreducible magnitude error Am: 



/0.15 2 

Nsn 



+ (Am) 2 



(10) 



As in Kim et al. (2004), we consider two different possibili- 
ties for Am: a constant value equal to 0.02 mag, which is the 
target of the SNAP mission, and a linear function increasing 
with redshift reaching at the maximum covered redshift, z = 
1.7, the maximum expected error for SNAP, 0.02 mag: Am = 
0.02(z/1.7). The resulting constraints on J7,„o and a are shown 
in Figure 7 for the RP potential in the case of OQ and EQ 
(top-left and bottom-left panels, respectively). In both cases 
we show the la confidence regions when no priors are con- 
sidered. The solid contour refers to the constraints obtained 
by considering statistical errors only; they are the same shown 
in the top-left panel of Figures 3 and 5. Dotted and dashed 
lines present the 1 a regions when we include constant and lin- 
ear irreducible systematic errors. We notice that systematics, 
as expected, extend the confidence regions, in particular in the 
direction of lower values for i7„,o and higher values for a. This 
effect is larger when we apply a constant Am. In the OQ case, 
the best-fitting values for a is shifted up of 0. 1 and the errors 
on the two variables are strongly correlated. 

Then we investigate the effect of systematic magnitude off- 
sets. First, we consider a constant shift of 0.02 mag, both pos- 



itive and negative. The analysis (not shown in the figure) con- 
firms the results obtained by Kim et al. (2004): in this case 
the best-fitting parameters are the true values with errors which 
are very similar to the ones obtained by considering statistical 
errors only, without any systematic biases. Then we consider a 
magnitude offset which is linearly proportional to the redshift 
as Am = ±0.03(z/1.7). The results for the same models previ- 
ously considered (OQ and EQ with RP potential) are shown in 
the right panels of Figure 7. As discussed by Kim et al. (2004), 
an offset in magnitude can give best-fitting parameters which 
are wrong with smaller confidence regions, i.e. it is possible 
to have very accurate but wrong answers. This is exactly what 
is happening in our case. Considering a positive linear offset, 
we find the tendency to have confidence regions systematically 
shifted towards the left, while we have the opposite trend for 
negative magnitude shifts: the true values of the parameters f2 m o 
and a are excluded at lcr confidence level. Finally we consider 
the presence of a possible mismatch between the calibration 
of low- and high-redshift SNIa. This is done by introducing 
a constant offset of ±0.02 mag to the 300 local objects only. 
The results (not shown in the figure for clarity) show that this 
kind of effect produces a small enlargement of the confidence 
regions, without introducing systematic biases in the parameter 
determination. 

5. CONCLUSIONS 

The main goal of this paper has been to discuss the possi- 
bility of using SNIa distances and redshifts to obtain reliable 
constraints on the parameters defining the quintessence mod- 
els. In particular we considered extended quintessence models 
with Ratra-Peebles potential, and, for completeness, ordinary 
quintessence models with both Ratra-Peebles and SUGRA po- 
tentials. 

As a first step, we studied the constraints which result from 
the analysis of the largest SNIa sample presently available (176 
objects, coming from Tonry et al. 2003, Knop et al. 2003 and 
Blakeslee et al. 2003). To this purpose, we used the exact red- 
shift evolution of the equation of state w(z), as numerically de- 
termined, avoiding any approximation. Our results show that 
for ordinary quintessence with a SUGRA potential it is not pos- 
sible to obtain significant limits on the potential exponent a, be- 
cause of the weak dependence of the equation of state on it. For 
ordinary quintessence with a Ratra-Peebles potential, it is pos- 
sible to obtain an upper limit on a: at the lcr confidence level 
we find a < 0.83. Our results, which are in agreement with 
a similar analysis made by Podariu & Ratra (2000) who used 
an older SNIa dataset, are consistent with a cosmological con- 
stant model, for which a = 0. We obtained a similar constraint 
for a when considering the extended quintessence models with 
Ratra-Peebles potential. Unfortunately, it is not possible to ob- 
tain useful constraints on the non-minimal coupling parameter 
£ between the scalar field and gravity. 

We then discussed the potential improvement on the previ- 
ous results when future SNIa samples, with a larger number of 
objects and more extended redshift coverage, will be available. 
To this purpose, we simulated SNIa datasets in different cosmo- 
logical models, with the characteristics of the expected SNAP 
satellite observations (almost 2000 objects up to z = 1.7). 

For ordinary quintessence with a SUGRA potential, we found 
that it will still be difficult to constrain the parameters even 
with the SNAP SNIa data. Considering models with the Ratra- 
Peebles potential, both in extended and ordinary quintessence, 




FIG. 7. — Effects of systematics errors on the lcr confidence regions for 
the parameters Q m0 and a for the ordinary quintessence model with a Ratra- 
Peebles potential (upper panels) and extended quintessence model with a 
Ratra-Peebles potential (lower panels). The results are obtained with the sim- 
ulated SNAP samples (datasets 1 and 3, respectively). Left panels show the 
effects of irreducible systematic uncertainties: dotted (dashed) contours refer 
to confidence regions obtained by including constant (linearly increasing) error 
Am (see text for details). Right panels present the effects of magnitude off- 
sets: dotted (dashed) contours refer to a positive (negative) shift of 0.02 mag. 
In all panels, the solid curve represents the constraints obtained by assuming 
the presence of intrinsic statistical errors, only. 



our results suggest that a can be determined with an error ^0.6 
(at the lcr significance level), while Q m o will be constrained 
with an error of approximately 0.03 (always at 1 a). In extended 
quintessence models, even by imposing the upper bounds on the 
time variation of the gravitational constant, it will be possible to 
obtain only an upper limit on the coupling constant: £ < 0.028. 

As a final issue, we discussed the systematic errors on the pa- 
rameter estimates which originate if a constant or linear approx- 
imation is used for the equation of state, w(z), instead of the ex- 
act redshift evolution. For all the considered quintessence mod- 
els, the confidence regions obtained with these approximations 
look narrower than the exact ones. As a consequence, these ap- 
proximations on w(z) lead to a systematic underestimate of the 
errors. Nevertheless, the set of cosmological parameters deter- 
mined by the fitting procedure are consistent with that used in 
the simulations. This means that the systematic errors induced 
by the assumed approximations are still smaller than the preci- 
sion allowed on the cosmological parameters by SNAP, even at 
the la confidence level. 
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